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Abstract 



Quantum phase estimation is one of the key algorithms in the field of quantum computing, but up 
until now, only approximate expressions have been derived for the probability of error. We revisit these 
derivations, and find that by ensuring symmetry in the error definitions, an exact formula can be found. 
This new approach may also have value in solving other related problems in quantum computing, where 
an expected error is calculated. Expressions for two special cases of the formula are also developed, in the 
limit as the number of qubits in the quantum computer approaches infinity and in the limit as the extra 
added qubits to improve reliability goes to infinity. It is found that this formula is useful in validating 
computer simulations of the phase estimation procedure and in avoiding the overestimation of the number 
of qubits required in order to achieve a given reliability. This formula thus brings improved precision in 
3 ■ the design of quantum computers. 



^ ■ Introduction 

>'• 

00 ■ Phase estimation is an integral part of Shor's algorithm [1] as well as many other quantum algorithms [2], 

designed to run on a quantum computer, and so an exact expression for the maximum probability of 
error is valuable, in order to efficiently achieve a predetermined accuracy. Suppose we wish to determine 
a phase angle <j) to an accuracy of s bits, which hence could be in error, with regard to the true value of 



4>, by up to 2 s , then due to the probabilistic nature of quantum computers, to achieve this we will need 



to add p extra qubits to the quantum register in order to succeed with a probability of 1 — e. Quantum 
registers behave like classical registers upon measurement, returning a one or a zero from each qubit. 
Previously, Cleve et al. [3] determined the following upper bound: 



PC = 



(i) 



Thus the more confident we wish to be (a small e), for the output to achieve a given precision s, the more 
qubits, p, will need to be added to the quantum register. Formulas of essentially the same functional form 
as Eq. ([1]), are produced by two other authors, in [5] and 0], due to the use of similar approximations 
in their derivation. For example, we have p = [log 2 ( + 2) + log 2 7r] , given in [3]. As we show in the 
following, these approximate error formulas are unsatisfactory in that they overestimate the number of 
qubits required in order to achieve a given reliability. 

The phase angle is defined as follows, given a unitary operator U, we produce the eigenvalue equation 
U\u) = c 2m< ^\u), for some eigenvector \u), and we seek to determine the phase (f> € [0, 1) using the quantum 
phase estimation procedure [5]. The first stage in phase estimation produces, in the measurement register 
with a t qubit basis {\k}}, the state [5] 

2*-l 

|0)sta g cl = ^E e27ri0fc | fc )- ( 2 ) 
k=0 
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If <fi = 6/2* for some integer 6 = 0, 1 , ... 2* — 1, then 

2*- 1 e 27ribfc/2 4 
l</>)stagcl = ' With Vk = 2 t/ 2 ' ( 3 ) 

is the discrete Fourier transform of the basis state |6), that is, the state with amplitudes Xk = Skb- We 
then read off the exact phase cf> = 6/2* from the inverse Fourier transform as |6) = J^\(j))- 

In general however, when <p cannot be written in an exact t bit binary expansion, the inverse Fourier 
transform in the final stage of the phase estimation procedure yields a state 

|0) = J"t|0) stagcl , (4) 

from which we only obtain an estimate for </>. That is, the coefficients Xk of the state |</>) in the t qubit 
basis {|fc)} will yield probabilities which peak at the values of k closest to <fi. 

Our goal now is to derive an upper bound which avoids the approximations used in the above formulas 
and hence obtain a precise result. 



Analysis 

In order to derive an improved accuracy formula for phase estimation, we initially follow the procedure 
given in [3], where it is noted, that because of the limited resolution provided by the quantum register 
of t qubits, the phase </> must be approximated by the fraction where 6 is an integer in the range 
to 2* — 1 such that 6/2* = 0.6i . . . 6 t is the best t bit approximation to 4>, which is less than (j). We then 
define 

5 = (f)-b/2 t , 

which is the difference between cf> and 6/2* and where clearly < <5 < 2~*. The first stage of the 
phase estimation procedure produces the state given by Eq. @. Applying the inverse quantum Fourier 
transform to this state produces 

2 4 -l 

10) = J>* \k) , (5) 

fe=0 

where 

1 2t_1 » 1 1 P 27ri2 t ,5 

e=o i — e 2 

Assuming the outcome of the final measurement is m, we can bound the probability of obtaining a 
value of m such that \m — b\ < e, where e is a positive integer characterizing our desired tolerance to 
error, where m and 6 are integers such that < m < 2* and < 6 < 2*. The probability of observing 
such an m is given by 

e 

pr{\m-b\<e)= £ \x b+e \ 2 . (7) 

i=-e 

This is simply the sum of the probabilities of the states within e of 6, where 

1 i-e 2 " 2ta 

x b+i — 2 i 1 _ e 27ri{<5-£/2*) ' ^> 

which is the standard result obtained from Eq. ([B]), in particular see equation 5.26 in [2], Typically at 
this point approximations are now made to simplify xi, however we proceed without approximations. We 
have 

2 _ 1 - cos(2^2*J) 

]Xb+el ~ 2 2 * l-cos(27r(<5-^/2*)) ■ [ ' 
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Suppose we wish to approximate </> to an accuracy of 2 s , that is, we choose e — 2* s 1 = 2 P 1 , using 
t = s+p, and if we denote the probability of failure 

e=p(|m-6| > e), (10) 

then we have 

l-cos27r2*<5 2 ^ 1 
£ ~ 2^ l-cos27r((5-£/2*)' ( ' 

This formula assumes that for a measurement m, we have a successful result if we measure a state either 
side of b within a distance of e, which is the conventional assumption. 

This definition of error however is asymmetric because there will be unequal numbers of states summed 
about the phase angle (f> to give the probability of a successful result, because an odd number of states is 
being summed. We now present a definition of the error which is symmetric about <j). 



Modified definition of error 

Given an actual angle <j) that we are seeking to approximate in the phase estimation procedure, a mea- 
surement is called successful if it lies within a certain tolerance e of the true value 4>. That is, for a 
measurement of state m out of a possible 2* states, the probability of failure will be 





m 


1 2tt\ 


( 










> 2 ~¥) 



(12) 



Thus we consider the angle to be successfully measured accurate to s bits, if the estimated <f> lies in the 
range (j>± ||f . Considering our previous definition Eq. ([TO]) , due to the fact that b is defined to be always 
less than </>, then compared to the previous definition of e, we lose the outermost state at the lower end 
of the summation in Eq. (fTT|) as shown in Fig. (p}. For example for p = 1, the upper bracket in Fig. (JlJ 
(representing the error bound) can only cover two states instead of three, and so the sum in Eq. (fTTj) will 
now sum from to 1, instead of —1 to 1, for this case. 
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P = 3 , E-3 



Figure 1. Defining the limits of summation for the phase estimation error. For the cases 
p = 1, 2, 3, we show the measurements which are accepted as lying within the required distance of </>, 
shown by the vertical arrow, which define the limits of summation used in Eq. (|13[) , 



Nielsen and Chuang [2] in the preliminary to Eq. 5.35, appear to have written incorrectly 2 P — 1 instead of 2 P 1 . 
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An optimal bound 

Based on this new definition then for all cases we need to add 1 to the lower end of the summation giving 



1-cos27t2«<E 2 ^ 1 (n] 

e 2 2t e= _i" 1+1 1 - cos M s ~ e / 2t ) 

and if we define a = 2 l 8 and rearrange the cosine term in the summation we find 

2^- 1 

1 — cos 2ira v-^ , n , , . 

t= l - 22t+ i E csc 2 -(a-^). (14) 

l=-2P- 1 + l 

Next, we demonstrate that the right hand side of Eq. (|14[) takes its maximum value at a = i. Since 
we know < a < 1, and since we expect the maximum value of e = e(a, t,p) to lie about midway between 
the two nearest states to generate the largest error, that is at a = 1/2, we will substitute a = | + A, 
where A <C \- To maximize e we need to minimize 

2 P — 1 

cos27r(i+Aj Yu csc2 J(V £ + A )' ( 15 ) 

V / l=-2P- 1 + l ^ ' 

as a function of A. Expanding to quadratic order with a Taylor series, we seek to minimize 

(1 - ir 2 A 2 + 0(A 4 )) (c + Cl A + c 2 A 2 + c 3 A 3 + 0(A 4 )) , (16) 

where Cj are the coefficients of the Taylor expansion of cosecant 2 in A. We find by the odd symmetry of 
the cotangent about L=\ that 

d = f E cotJ(I-^)csc 2 J(i-£) = 0, (17) 

<=-2*- 1 +l 

and so we just need to minimize 

c Q + (c 2 -c ir 2 )A 2 + O(A 3 ). (18) 

Differentiating, we see we have an extremum at A = 0, and therefore e(a, t,p) has a maximum at a = 1/2. 
Substituting a = \ we obtain 

£ i-cos^-O - (19) 

We note that the summation is symmetrical about I = 1/2, and substituting t — p + s, we obtain for our 
final result 

op— 1 

1 1 

^P) = 1 - ^PR E 1 ( 2 °) 

That is, given a desired accuracy of s bits, then if we add p more bits, we have a probability of success 
given by 1 — e, of obtaining a measurement to at least s bits of accuracy. Thus we have succeeded in 
deriving a best possible bound for the failure rate e = e(s,p). 
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Special Cases 

Numerical calculations show that e(t,p) quickly approaches its asymptotic value as t — > oo, and this limit 

2 

gives a fairly accurate upper bound for e, for t greater than about 10 qubits. Using cosx > 1— %■ which 
is valid for all x, and is accurate for x = 0(1/2*) as t — > oo, 



e 1 - 



9 p — 1 

4 



y L 



2 2 t l_(l_|(-(2^-l))2) 



2 P— 1 

= '-^g^v (21) 

An exact form for this can be found in terms of the trigamma function, being a special case of the 
polygamma function as shown in Abramowitz and Stegun [6], Eq. 6.4.5 : 

where ip f (z) = ^7 is the trigamma function, ip(z) = YTz) * s ^ e digamma function, and T(z) = 
Jo°° t z ~ le t dt is the standard gamma function. 

Now considering the p — > 00 limit, which also includes the t — > 00 limit because t = p + s, we can find 
an asymptotic form in the limit of large p also from [5], Eq. 6.4.12, namely 

e = 4^, (23) 

which shows that the error rate drops off exponentially with p extra qubits. The formula Eq. (|23|) can be 
re-arranged to give 

2V2 

log 2 ■ 

which can be compared with the previous approximate formula shown in Eq. ([1]). 



(24) 



We have checked the new error formula through simulations, by running the phase estimation algo- 
rithm on a 2-dimcnsional rotation matrix, and undertaking a numerical search for the rotation angle that 
maximizes the error e, which has confirmed Eq. (|20[) to six decimal places. 



Discussion 

An exact formula is derived for the probability of error in the quantum phase estimation procedure, as 
shown in Eq. (|20[) . That is, to calculate <p accurate to a required s bits with a given probability of success 
1 — e we add p extra qubits, where p is given by Eq. (|20[) . If we have a large number of qubits then 
we can use Eq. (|22[) valid at the t — > 00 limit. In the p —> 00 limit the asymptote is found as a simple 
exponential form Eq. (|2"5]l . 

The exact formula avoids overestimating the number of qubits actually required in order to achieve a 
given reliability for phase estimation and we have also found this formula to be useful in confirming the 
operation of classical simulators of the phase estimation procedure. 
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